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ABSTRACT 


Thermoelectric Cooler Design 1s a Microsoft Windows program to be used as an aid in 
the design of thermoelectric cooler devices. This program was written to be used to 
quickly model and compare alternative designs. A couple's optimum coefficient of 
performance and maximum heat pumping can be quickly determined. Other major 
features of the program include the ability to change material properties and dimensions 
of couples, analyze cascaded couples, and graph performance parameters. A brief 
description of thermoelectric cooler theory, modelling assumptions and complete source 


code listing is included. 
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I. INTRODUCTION 


A. FEATURES OF THERMOELECTRIC COOLERS 


The unique features of thermoelectric coolers include: 


l. Compact. 

No moving parts. 

The coefficient of performance is independent of capacity. 

The cooling rate is easily controlled. 

Reversible operation is possible allowing heating as well as cooling. 


Vibration-free operation. 


P ኮረ A rur 


Precise temperature control. 


B. SPACEBORNE APPLICATIONS 


Thermoelectric coolers offer an alternative solution to cooling and heating 


requirements of spacecraft. Spaceborne applications for thermoelectric coolers include: 


l. Equipment cooling -- removal of waste heat in order to maintain 
components at their optimal operating temperature. Thermoelectric coolers 


are especially useful for the control of small hot spots. 


2. Life support temperature control -- thermoelectric coolers are very reliable 
making them ideally suited for long duration missions such as the space 
station or manned Mars mission. 

3. Cooling of infrared detectors -- thermoelectric coolers eliminate the need 
for coolant and associated support equipment. Of even more importance is 
the elimination of consumable coolants which also increases the life of the 
spacecraft. 

4. Microgravity applications -- thermoelectric coolers offer precise, 
vibration-free temperature control for free-flying orbiting microgravity 


research platforms. 


Reliability, the lack of moving parts or consumables, and vibration-free operation 


are especially important to spaceborne cooling systems. 


Il. THERMOELECTRIC COOLER THEORY 


A. THERMOELECTRIC COUPLE CONFIGURATION 


In it's most basic form, a thermoelectric couple consists of two dissimilar elements 
with contact at two junctions. Figure 1 shows a simple thermoelectric circuit consisting 


of two elements, physically joined on one side while electrically joined at the other. 





Figure 1 Thermoelectric Cooler 


B. THERMOELECTRIC EFFECTS 


The background for this section was gathered from Kraus [Ref. 12:pp. 294-303] and 
Ioffe [Ref. 2:pp. 96-128]. Three thermoelectric effects are known; the Seebeck, Peltier, 


and Thomson effects. 
1. Seebeck Effect 


Consider the junction of two materials, A and B. If a temperature differential 
exists between the two junctions, one Junction may be referred to as the hot side and the 
other the cold side. An open circuit voltage can be measured between the two junctions. 
This net conversion of thermal energy into electrical energy under zero current 
conditions is known as the Seebeck effect. The voltage observed, referred to as the 
Seebeck voltage, is proportional to the temperature differential between the hot and cold 


sides of the couple: 
dEs = dT 


dE ved 


Es = 1 | adr (V) (1) 


where 


Es = Seebeck voltage (V) 
ao = Seebeck coefficient of material (V /^C) 


In the case of the couple composed of materials A and B, 


œ= [oy | ር [ag 


where 


Ay = Seebeck coefficient of material A 
ap = Seebeck coefficient of material B 


2. Peltier Effect 


The Peltier effect is localized at the junctions between materials 4 and B. When 
a current flows through the junction, heat is either generated or absorbed depending on 
the direction of the current flow. The amount of heat is proportional to the current and is 


known as the Peltier heat: 
d Qp ec [dt 


d Qp —- tmnldt 


gp = intl (W) (2) 
where 


Qp - Peltier heat (J) 

qp = Peltier heat flow (W) 

I = Current (A) 

n = Peltier coefficient of material (V) 


3. Thomson Effect 
This effect concerns the absorption or generation of heat in a conductor carrying 
a current in the presence of a temperature gradient. The amount of heat is proportional to 


the current and temperature gradient: 
dO, = IdTdt 


dO, = tol dTadt 


dq, SEO 


T) 
q; 7 tol|dT (W) (3) 
7] 


where 


Q; = Thomson heat (J) 
gt = Thomson heat (W) 


I = Current (A) 
o = Thomson coefficient of material (V/°C) 


Because the product of o and AT is a voltage, the Thomson voltage is therefore: 


T) 
EN T (4) 
T 


C. IRREVERSIBLE EFFECTS 


1. Joule Effect 


The evolution of heat when a current passes through a material may be referred 


to as the Joule heat flow: 


qj =R (W) (5) 


In the case of a couple, the elements 4 and B are connected in electrical series 


regardless of the direction of the current flow: 


where 
P4>PB 
A 
A4, AB 


L T 
Á 4 AB 





electrical resistivity of materials Æ and B. (€) cm) 
Length of materials A and B (cm) 


= Cross sectional area of materials A and B (cm?) 


2. Fourier Effect 


The conduction of heat through a material is given by: 


a= (KA Jar cw (7) 


and, in the case of a couple with the elements in thermal parallel, 


where 


k 4, kp - 
L4; Lg 
pun 





kada  kpA 
K= ሥማ ታር.) (8) 
L4 Lg 


Thermal conductivity of materials Æ and 8. (W/cm ^C) 
Length of materials A and B (cm) 


Cross sectional area of materials A and B (cm? ) 


D. NET HEAT ABSORBED AT COLD JUNCTION OF COUPLE 


In this analysis, the Thomson heat will be neglected which greatly simplifies couple 
design procedure without any significant error Kraus [Ref. 1:p. 300] and Cadoff [Ref. 
3:p. 8]. To determine the net heat pumped by the couple, the temperature distribution 


along the elements of the couple must be determined. Consider a single thermoelectric 
element as indicated in Figure 2, where Ax is the cross sectional area of the element 


normal to the flow of current and heat. 





Figure 2 Thermoelectric Element 


Here, the boundary conditions are: 


TS To at x= 
l= T2 txFL 


Assuming steady state temperature conditions, there are three heat quantities which 
pertain to an energy balance over the incremental length, Ax: 


1. The Fourier heat entering at point x: 


e ar 
ተ 6491). 


2. The Joulean heat generated in the increment Ax : 


B I? p Ax 
m Áx 





3. The Fourier heat leaving at point x * Ax : 


dT 


> 0 $e 
1 x+Ax "| dx x+Ax 


The energy balance in the steady state requires that the heat leaving must be equal to 


the sum of the heat entering and the heat generated in the increment Ax : 
dx * dj 7 d xe Ax 


Substituting for gx, 9 y+Ay and qj gives 


Dividing through by Ax yields 


EMT 
dx Ix+Ax dx JS x I^ p 





k Ax | A tt n 0 
Ax 
And taking the limit as Ax — 0 gives 
im, (len ER): 
1mit 
e0 A Al NE 
Áx 
Ax 
which can be adjusted to 
2 
d^T | Ip. 
kAx ው + A = 0 (9) 


The general solution to this simple differential equation is obtained through the use of 


a double integration: 





2 2 
Er i መ+ር1=0 
dx kAt 
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and 
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E Po + Cyx+Cy =0 
2k AZ 





Applying the boundary conditions for the hot side where x = O0 and T - To: 


"nM 


and 


I? px 


2kA2 


2 
T — To + 





+Cıx=0 


At the cold side of the element where x = L and T = Te: 


2 


I pL? 
Te — To + x +CıL=0 
2kA* 





The solution for Cy, the first constant of integration, 1s seen to be 


lo- lem ion 


Cy = 
L 2242 


so that the particular solution of this simple, second order, differential equation which 


gives the temperature profile as a function of distance along the element: 








2 L 2 


29 i 2 
um ES » (ferie). - |— k= 0 
2kAx 2ጽ4ሄ 


Or 








2 2 
I _ 
Eee — x2 - | Tote mp can (10) 
2k A 2k A 


which is seen to be parabolic. 


Figure 3 illustrates a temperature distribution for a typical element with the following 


parameters: 
231 MM 
z = 3.08x 10 (C ) (11) 
R2173x10^ (Q) (12) 
K = 3.64x1072 (W/°C) (13) 
ao = 44x10 (V/°C) (14) 
Tc 2 3086 K (15) 
AS IC (16) 


Due to the parabolic temperature distribution, the Joulean heat will be transferred to both 


the cold and hot junctions of the couple. 


Differentiating equation (10), setting the derivative equal to zero, 





ar_ BpX To-Te . pk 
dx ^ A ካዓ 
kA£ E. 


permits a solution for x 2 X which is the position where the maximum temperature 
occurs: 
FpL_To-Te 
244: d 
[2 pX 
ZE 


= 
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Figure 3 Temperature Profile Along Element 


Rearranging and substituting AT for To — Te: 


k AT AZ 


LI —— {| 


Assuming x = L and defining f as the fraction of the Joulean heat transferred to the 


cold junction: 
/- ፊ2ዳ Q1 4 KATAR 
Or 
28 ስአ 


where K and R are defined by equations (6) and (8), respectively. 
The net heat absorbed at the cold junction of the couple is equal to the summation of 


the Peltier heat, minus the fraction of Joulean heat transferred to the cold junction: 
One alter 
and upon substitution ofthe expression for f given by equation (17): 


gner = @ITec - SPR - KAT (W) (18) 
Therefore the net heat pumped is the Peltier heat developed at the cold Junction less 


the sum of one half the Joulean heat produced in the couples and the heat due to 


conduction from the hot junction to the cold junction. 
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III. THERMOELECTRIC COOLER MODEL 


A. INTRODUCTION 


The equations in this chapter were used to model the thermoelectric cooler couple 


operating parameters and were given by Kraus [Ref. 1:pp. 303 - 326]. 


B. HEAT PUMPING 


1. Current Yielding Maximum Heat Pumped 


The equation for net heat pumped 


aa Ao ¿PR - KAT (W) (18) 


may be differentiated with respect to current, /. When the derivative 1s set equal to zero, 


the current that yields maximum heat pumped is obtained: 


dq 


paa o. 0 
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Solving for I: 


Im = on (A) (19) 





2. Maximum Heat Pumped 


Substituting /m into the equation for the net heat pumped yields the maximum 
heat pumped: 
02 e 
2R 





qm = AI) (20) 


C. MAXIMUM TEMPERATURE DIFFERENTIAL 


Heat pumped as a function of current 15 shown in Figure 4 using the parameters 


defined in equations (11) through (16). Setting equation (20) which is the maximum heat 
pumped equal to zero and solving for AT, one obtains the maximum temperature 


differential between the hot and cold sides of the couple: 


2 ተረ 


OI 
Alm = Tan (K) (21) 


D. OPTIMUM RATIO OF ELEMENT CROSS SECTIONAL AREAS 


The maximum temperature differential given by equation (21), can be ıncreased 


by minimizing the product KR from equations (6) and (8): 
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aídno>/(m) b 


m 


ከዚህያም 
KEF 
— COS 


of Current 


Figure 4 Heat Pumped as a Function 


iL L k4 A kp A 
ka = (2474 + Pe le re 22) 
A4 AB La Lg 


andif L, = Lp an expansion gives 
AB 44 
KR O04 4 1 EDIP 7 (23) 


An optimization of equation (23) can be accomplished by finding the point 


where the derivative with respect to the area ratio, 4 4/4 p, vanishes: 


2 
. GKR — — => AB = 
d(A4lApg) - “APB ~ KBP (32) e (eu) 


Then, a solution of equation (24) yields the area ratio that minimizes KR. Thus, in turn 
yields the optimal AT m: 
A4 _ [ቦ4ጃጩ 
AB VPBKA 





(25) 


The area ratio given by equation (25) when substituted into equation (22), yıelds 


the optimum value of KR, designated by $ 


9—k4p4 *2/k4pA4kpppg * kpgpg 


Or 


+= (¡rana + [hapa ) (26) 
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Equation (16) has the restrictions that L 4 = Lp and that the area ratio be determined 


using equation (25). 


E. FIGURE OF MERIT OF ELEMENT MATERIALS 


Substituting $ , equation (26),into the equation for ATm equation (21) : 


2 
ር... 
OT 
ATm - LiT 


where z, the figure of merit of the materials, 1s defined as: 


T3452 
Q 


which may be rewritten as: 


a2 


(apa + ffo ) 


Z = 


When L 4 = Lp and the area ratio determined from equation (25), the figure of 


merit 1s a function of the properties of the elements. To maximize the figure of merit, the 
materials used should have a high Seebeck coefficient, low thermal conductivity, and low 


electrical resistivity Cadoff [Ref. 3: p. 21]. 


F. POWER REQUIRED TO OPERATE THE COUPLE 


To operate the couple, a voltage equal to the sum of the Seebeck voltages in the 
couple and the resistive voltage drop must be applied Kraus [Ref. 1: p. 304]. The power 


required to operate the couple is then: 


P = IV = IQ AT + IR) (W) (27) 


G. COEFFICIENT OF PERFORMANCE 


The coefficient of performance,n , measures the performance of the couple and is the 


heat pumped, g, divided by the power required to operate the couple, P: 


m. 
NT (28) 


Substituting equation (18) for q and equation (27) for P gives 


al Te - 51? R-KAT 
‘lee AIR + QAT) (29) 


1. Coefficient Of Performance At Maximum Heat Pumping 


One may express qm in terms of AT and ATm using equations (20) and (21): 
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ጄም 2 
Ka? T2 
ር መድ ለ 


Tc 
[ሀመ 7227 ርን 2 


qm = KATm - KAT = K(ATm = AT) 





qm = K ATm [1 x - | 


(30) 


In similar fashion, Pm may be expressed in terms of the temperatures by 


substituting equation (19) into equation (27): 
Pm = & /ImAT + IŻ R 


. GiTgaT r 


B R R 





. 2Ka? Té (1 AT) 


= — 2KR p 


RN. n 7 Ar) 
Te 


To obtain the coefficient 


(31) 


of performance when the couple is adjusted for 


maximum heat pumping, substitute equations (30) and (31) into equation (28) for the 


coefficient of performance: 


2] 


E qm å KATm (1 — AT/ATm) 
nm E INTA CENA 


AA 
Am = IT + ATIATE) 


2. Current Yielding Optimum Coefficient Of Performance 


The expression for the coefficient of performance, equation (29), 


a 1 Tc - 3I R-KAT 


I 


may be differentiated to find the current that yields the maximum coefficient of 


performance. This current is the value of Z that causes dn/dI to vanish: 


dy _ IUR+0AT) (QTc -IR)- (ITee -I R/2- KAT QIR +4AT) _ 


= 0 
dI I? (IR - . AT)2 
Or 
(PR * aJAT) (a Tc — IR)» (a1Tc — I R/2 — KAT) QIR + a AT) 
Expansion gives: 
PRaTe - DR? & a2ITc AT - aIRAT 2 2 Ra Te * 121 Tc AT 
-PR2 - au PRATI2 - 2IRKAT - AKAT? 
OT 


al2RAT/2 + a Tc I?R - 2KIRAT — KaAT* = 0 
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which is actually a quadratic in J: 


Ro (47 E Te) P EO RRN- (Kaar? ) = 0 


which possesses a solution 


2KRATt |Q KRAT)? - 4KRO?AT?[(AT/2) + Tc] 


Oooo GRaQTEmAWRb እእ ES 


The value of the radical in equation (32) will always be greater than 2KR AT 
because it contains the sum of (2፳፳ላ7)* and another, positive term. For cooling to 


occur, the current must not be negative so the solution in which the radical term is 


subtracted is ignored: 


2KRAT+ J(2KRAT)? + 4KRa2AT?[(AT/2) + Te] 


factoring 4 (KR)? AT? from under the radical glves: 


; 2KRAT + 2KRAT/1+ (02/KR)[(AT/2) + Tc] 


7 2 Ra [(AT/2)+ Tc] 


or 


1 + /1 + 2 [(ላ7/2) + Tc] 


ከጨ Ra [(AT/2) 4 Tc] 
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Because 


DE pe 
AT t Te = A + Te = 4 





we define the average temperature may be defined as: 





_ Ty+T 
ys ው - (33) 
and 
po Karp Aita *TO?2] - 
E (Ty + ፲ር)/2 (34) 


Multiplying equation (34) by 


iue zd 
RO! Ji+zT -1 


and substitution of the average temperature given by equation (33) one obtains: 


I 


E ERST 4 qure 
Ro? T | ገ TE] 


With 1/2 = KR/o?, expansion gives: 


T ser zT | 
R ze ea 
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Rearranging yields the expression for the current that optimizes the coefficient of 


performance: 
፪ ላያ 


ጸ[ሀ1+2፻ -!] 


lo = (35) 


Figure $ 1s a plot of the coefficient of performance as a function of current based 
on the conditions specified in equations (11) through (16). 


3. Optimum Coefficient Of Performance 


Substituting Jo into the coefficient of performance equation (35) yields the 


optimum coefficient of performance: 


ገ በ ም Ih 
ቸም ር. ጠሙ ክሮ 
AT Jl+zT +1 


(36) 


As z approaches œ, No approaches the ideal, or Carnot, coefficient of 


performance which is the term to the left of the brackets in equation (36) Egli [Ref. 4: p. 


31] and Tipler [Ref. 5:pp 575-576]. 


H. CASCADED COOLERS 


If a large AT is required, a multistage or cascaded cooler can be used. In a cascaded 


cooler the heat rejected by the first stage of couples is the heat load fed to the second 
stage. This arrangement can be extended to any number of stages. For best results when 


designing a multistage cooler for the lowest temperature attainable, each succeeding 
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Figure 5 Coefficient of Performance as a Function of Current 


stage should require less power than the previous stage. It is also possible to increase the 
coefficient of performance when the required temperature drop is close to the maximum 


for a given figure of merit Ioffe [Ref. 2: p. 116]. 
1. Overall Heat Pumped By Cascaded Coolers 


The heat pumped by a battery of n-stages of cascaded coolers 1s equal to the 


summation of the heat load, q at the first stage and the power required, P, of each stage 


over all stages in the battery: 


N 
doverall 7 41* 2 Pi (W) 


i= 


2. Overall Coefficient Of Performance 


The coefficient of performance of a battery n-stages of cascaded coolers is the 


heat load of the first stage, q, , divided by the total power required to operate the battery 


of coolers. 
d] 


N 
ye Wee 
isl 





Noverall = 


I. SUMMARY OF SIMPLIFYING ASSUMPTIONS 


1. Thomson Voltage 


Neglect of the Thomson voltage results in no significant error. However, designs 
utilizing this model will be conservative, understating couple performance Kraus [Ref. 


| 0 1 
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2. Temperature Losses 


Temperature losses between the thermoelectric couple and the heat load or heat 
exchanger are not taken into account directly. However, once determined, these 
temperature losses can be compensated for by adding the drop to the cold side 


temperature and subtracting the drop from the hot side temperature. 
3. Parallel Thermal Paths 


Insulation between couples is assumed to be perfect, eliminating parallel thermal 


path heat losses. 
4. Junction Resistance 
Junction resistance is assumed to be zero. 
S. Figure Of Merit 


In order to maximize the figure of merit, the couple material lengths are assumed 


to be equal and area ratio determined by equation (25). 


6. Resistivity And Conductivity 


Electrical resistivity, p , and thermal conductivity , k, are assumed to be constant 


with respect to temperature. 
7. Couple Power Supply 


The couple power supply is assumed to be direct current. 
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IV. THERMOELECTRIC COOLER DESIGN PROGRAM 


A. INTRODUCTION 


1. Features 


Thermoelectric Cooler Design is a Microsoft Windows program to be used as an 
aid ın the design of thermoelectric cooler devices. Use of the Microsoft Windows 
graphical user interface allows the user who is familiar with the Windows environment to 
concentrate on the design process instead of a new, nonstandard interface. This program 
was written to be used to quickly model and compare alternative designs. Modifications 


to a design can be easily incorporated. Major features of the program include: 


Ability to change material properties and dimensions of couples. 
Determine optimum coefficient of performance of each couple. 
Determine maximum heat pumping capacity of each couple. 
The ability to analyze cascaded couples. 


Graph the performance parameters of couple. 


Nn (ቂል A UU N -e 


Provide on-line help in the form of a glossary, diagrams, and a temperature 


conversion calculator. 
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2. System Requirements 


The minimum software and hardware requirements needed by the computer 


system to run Thermoelectric Cooler Design are: 


l. Microsoft Windows, version 3.0 or later. 
2. A hard drive with 500 kilobytes of free space. 


3. A mouse that is supported by Windows. 


3. Disk Contents 


Files included on the installation disk: 


l. VBRUNIOO.DLL - Visual Basic kernel which must be in the same 
directory as THERMAL.EXE. 

2. THERMAL.EXE - Thermolelectric Cooler Design executable file. 

3. README.TXT - This manual in ASCII format. 


Both VBRUNIOO.DLL and THERMAL.EXE must be installed in order to run 


Thermoelectric Cooler Design. 
4. Installation 


To install Thermoelectric Cooler Design: 


l. Insert the installation disk in A: drive. 

Exit to DOS. (Ensure C:\ prompt 1s present). 

Type MD THERMAL then press the ENTER key. 

Type COPY A:\*.* C:\THERMAL then press the enter key. These 


E E 


commands create a new directory named THERMAL on the C drive and 
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l. 


copy the files from the installation disk into the C:\THERMAL directory. 
However, this is only one of many different ways to copy the files from the 
installation disk to the computer. If the choice is to install Thermoelectric 
Cooler Design on another drive or in another directory, both 
VBRUNIOO.DLL and THERMAL.EXE must be installed together in the 
same directory. The program will not run if this is not the case. 

Run Microsoft Windows. 

Select the Program Manager. 

Open the Group Window in which Thermoelectric Cooler Design is to be 
installed. 

Choose the option "new" from the file menu of the Windows Program 
Manager window. The new program object dialog box will then appear. 

Select the program item and choose OK. The program item properties 


dialog box appears. 


. Type a description to appear under the program icon. 


. Select the browse button to display a list of files and directories. 


Select THERMAL.EXE and choose OK. 
Select OK to close the program item properties dialog box. 


The Thermoelectric Cooler Design program is now installed. 


B. USING THE PROGRAM 


Starting And Exiting The Program 


To run Thermoelectric Cooler Design, simply double-click the program icon. 
After displaying the introduction, the program menu bar will appear. This program menu 
bar, titled "Thermoelectric Cooler Design" provides access to all program functions. To 
increase readability, select the minimize on use option under the Windows program 


manager options menu. 


3] 





Thermoelectác Cooler Design it. 
. Glossary Diagrams Material Properties | Optimize Graph Calculator Exit 


Info 
Figure 6 Program Menu Bar 


To exit the program, select exit from the program menu bar and choose yes 


when prompted. Once exited, all current design data ıs lost. 
2. Setting Material Properties 


Prior to determining the optimum coefficient of performance or maximum heat 
pumping capacity of a couple, the material properties and dimensions of the couple must 
be input. 

To begin, select material properties from the program menu bar. The material 
properties menu will appear with default values. If these values are acceptable, choose 
OK and proceed to the optimize menu. If the default values are not acceptable they may 


be changed by: 


1. Choosing either the first or second material. The «, p and k of the selected 
material is then displayed. All values in blue can be changed by the user. 
Appropriate values of &, p and k may be entered by the user by highlighting 


the text to be changed with the mouse and inputting the desired value from 
the keyboard. 

2. Changes to the dimensions of the couple are input in a similar manner. Both 
the area and length of the second material may be input directly by the user. 

3. Once the desired material properties and couple dimensions have been 
input, the overall resistance, thermal conductivity, and figure of merit are 
calculated by selecting the update button. In order to proceed to the optimize 
menu, select OK. To return all values to the defaults, choose the default 


button. 
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Properties Dimensions 
@ First Material 


Second Material 
- Material 
2 
œ 10.00023 | v79 c Årea 0.3554 | cm 


Effective 
D 10.00100! Q cm Diameter 0.6727 | cm 0.7000 | cm 


k a n cm cm 


Overall 


Resistance (A) Q 
Thermal Conductivity (K) [3.659E-02_|W/ °c 


Figure of Merit (2) 


Figure 7 Material Properties Window 





3. Optimum Coefficient Of Performance 


To calculate the optimum coefficient of performance of a couple: 


1. Ensure that the material properties menu has been reviewed. 
2. Select the coefficient of performance option of the optimize menu. 


3. The performance window will then appear. For a given low temperature the 
maximum high temperature is computed. The default AT ıs the maxımum 
and may be changed by manipulating the vertical scroll bar. Once the heat 
load, low temperature and AT are acceptable, the optimum coefficient of 
performance is calculated by pressing the compute button. The optimum 


coefficient of performance and associated performance values are then 


displayed. Changes to the heat load, low temperature, and AT can be made 


and all couple performance values recomputed by choosing restart. 
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Figure 8 Optimum Coefficient of Performance Window 
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Figure 9 Coefficient of Performance Results Window 
4. Maximum Heat Pumping 


To calculate the maximum heat pumping capacity of a couple: 


1. Ensure that the material properties menu has been reviewed. 
2. Select maximum the heat pumping option of the optimize menu. 
3. The pumping window will appear. For a given low temperature the 


maximum high temperature is computed and displayed. The default 
temperature difference, AT, is the maximum and may be changed by 
manipulating the vertical scroll bar. Once the heat load, low temperature, and 
AT are acceptable, the maximum heat pumped is calculated by choosing 
compute. The maximum heat pumped and associated couple performance 
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values are then displayed. Changes to the heat load, low temperature and AT 
can be made and all couple performance values recomputed by choosing 









restart. 
Load (W) 
4552 
High 
Temp {K} Delta T: 
455.1 146.5 
Low 
Temp {K) 
308.6 308.6 
Figure 10 Maximum Heat Pumped Window 
Print 
A/couple 300.0 | Number of Couples 


W /couple Coefficient of Performance |g 3294 


V/couple [9.1525 | 10.7253 Total Heat Rejected 12158 |W 
15.3626 | 
11.3941 _ Power C ü W 
ower Consumption 


Figure Eleven Heat Pumped Results Window 


5. Cascading Stages 


Once the initial couple performance values have been computed for either 
optimum coefficient of performance or maximum heat pumping, cascading of multiple 


couple stages can be simulated by selecting cascade. The cascade button is not visible 
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until the first stage couple performance has been calculated. Each additional stage 
assumes the rejected heat of the previous stage at the heat load and the high temperature 


of the previous stage. The maximum temperature at the new low temperature is then 


calculated. While the AT of the previous stage is assumed as a default, the AT for every 
stage of the cascade can be changed with the vertical scroll bar. While the program only 
displays five stages at a time, an unlimited number of cascade stages can be computed 
and displayed. Results of stages beyond five are shown in either the optimum coefficient 


of performance or maximum heat pumping windows, but not in the cascade window. 
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Figure Twelve Cascaded Stages Window 
6. Graphing The Coefficient Of Performance And Heat Pumped 


Thermoelectric Cooler Design is able to graph both coefficient of performance 


and heat pumped as a function of current. These graphs are useful in determining the 
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effect of a small change in current on both the maximum coefficient of performance and 
the heat pumped. The effects of varying the temperature differential across the couple 
can also be compared. 

To graph either the coefficient of performance or heat pumped as a function of 


current: 


1. Ensure that the material properties have been reviewed and optimum 
coefficient of performance or maximum heat pumped has been determined. 
The graph option of the program menu bar will not be available until these 
steps are accomplished. 

2. Select the graph option from the program menu bar. The graph window will 
appear. 

3. Select the function option from the graph window menu bar. Choose either 
coefficient of performance versus current or heat pumped versus current. 


4. Choose draw from the graph window menu bar. The selected function is 
then graphed for the temperature differential, AT, chosen in the performance 
or heat pumped window. 

S. To change the temperature differential, manipulate the AT vertical scroll 
bar in either the performance or heat pumped window as appropriate. 


To superimpose plots of representative of various temperature differentials and 
functions, select superimpose from under the options choice of the graph window menu 


bar. For example, to compare the effect of varying the temperature differential on a 
graph, change the AT vertical scroll bar in either the performance or heat pumped 


window. Subsequent draw commands are then superimposed on the current graph. To 


clear the plot, deselect superimpose and then select draw 
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Figure Thirteen Optimum Coefficient of Performance Graph 
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Figure Fourteen Maximum Heat Pumped Graph 
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The number of samples used to produce the graph may be adjusted as required. 


Select the sample size desired under the samples choice under the options selection of the 


graph 


Ot A oU 


window menu bar. 
7. Temperature Conversion Calculator 


To convert between temperature scales: 


Select the calculator option of the program menu bar. The temperature 
conversion calculator will appear. 

Input the temperature to convert from by highlighting the blue numerical 
display with the mouse by holding down the left mouse button while 
sweeping the cursor from left to right across the calculator display text. Once 
the text is highlighted, type in the desired temperature. 

Select the temperature scale to convert "from". 

Select the temperature scale to convert "to". 

Choose convert and read the temperature in the calculator display. 


To clear the display, select clear display option. 
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Figure Fifteen Temperature Conversion Calculator 
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8. Diagrams And Glossary 


On-line assistance is available in the form of diagrams and a glossary. 
Diagrams of a single couple and cascaded couples are available by selecting the 
diagram option from the program menu bar. 
A glossary of terms used throughout the program is available. Choose the glossary 


option of the program menu bar. 
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V. CONCLUSIONS 


The biggest drawback of a thermoelectric cooler is it's inherent low coefficient of 
performance, but new materials with better thermoelectric properties could eventually 
make thermoelectric devices competitive with compressor systems. While unable to 
compete directly with compressor based systems, thermoelectric coolers do have many 
unique characteristics making them useful in a variety of applications. For spaceborne 
applications, high reliability, vibration free operation, and no dependence on 
consumables are their most important features. 

The design of thermoelectric coolers 1s tedious without the aid of a computer. 
Unfortunately, the utility of many programs written for DOS are handicapped by a poor 
user interface. More time is often required to figure out how to use the program than to 
solve the initial problem. In writing thermoelectric cooler design, the Windows graphical 
user interface was used to provide a better user interface. Windows offers the 
programmer several powerful program development tools. For the user, Windows offers 


a familiar working environment enabling the user to concentrate on the task at hand. 
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APPENDIX 


Computer code listing for Thermoelectric Cooler Design program. 


Thesis.txt 
// Code for main menu bar 
Sub ProgramInfoName, Click () 
MsgBox "LCDR Clifton's Naval Postgraduate School Thesis" 
End Sub 


Sub DefinitionsName . Click () 
Definitions.Show 
End Sub 


Sub TemperatureName, Click () 
Screen. MousePointer = 11 
TempCalc.Show 
TempCalc. WindowState = 0 
Screen.MousePointer = 0 

End Sub 


Sub ChooseName Click () 
Material.Show 
End Sub 


Sub ShowGraphName Click () 
Screen.MousePointer = 11 
Graph.Show 
Graph.WindowState = 0 
Screen.MousePointer = 0 

End Sub 


Sub ExitName_Click () 
Screen.MousePointer = 1] 
Ending.Show 
Screen.MousePointer = 0 

End Sub 


Sub Form_Load () 
Main.Left = (Screen. Width - Main. Width) / 2 
Main.Top = 0 

End Sub 
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Sub Form Click () 
Main.Left ^ (Screen. Width - Main. Width) / 2 
Main.Top = 0 


End Sub 


Sub PerformanceName Click () 
End Sub 


Sub PumpingName Click () 
Pumping.Show 
End Sub 


Sub Commandl Click () 
Results.CoeffText.Text — "" 
Results. NumberCouplesText.Text = 
Results.TotalHeatRejectedText.Text — "" 
Results. AperCoupleText.Text — "" 
Results. VperCoupleText.Text = 
Results. WperCoupleText.Text = 
Results.qTText.Text = "" 
Results.qPText.Text ="" 
Results.aJText.Text="" 
Results.qFText.Text = "" 
Results.qoText.Text = "" 


"tt 


አዚ.) 


mo 


End Sub 
Sub InfoName_Click () 
MsgBox "LCDR Clifton's Naval Postgraduate School Thesis" 
End Sub 
Sub SinglePumpingName_Click () 
Screen.MousePointer = 11 
ScrollUpdate = False 
Te = Val(Pumping.LowTempText.Text) 
DeltaTMax 7» .5 *z * (Tc^2) 
TMax - DeltaT Max + Tc 


Pumping.Tempscroll.Max = Tc * 10 + 1 
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Pumping.Tempscroll.Min » TMax * 10- 1 
Pumping.Tempscroll.Value = Pumping.Tempscroll.Min 


Pumping.TMinDisplay.Text = Format$(Tc, "##,##0.0") 
Pumping. TMaxDisplay.Text = Format$(DeltaTMax + Tc, "##,##0.0") 


Th = Pumping. Tempscroll. Value / 10 
Pumping.HighTempDisplay.Text = Format$(Th, "##,##0.0") 
DeltaT = Th - Tc 

Pumping.DeltaTDisplay.Text = Format$(DeltaT, "##,##0.0") 


Pumping.Show 
Pumping. WindowsState = 0 


Screen.MousePointer = 0 
ScrollUpdate = True 
End Sub 


Sub GlossaryName_Click () 
Screen.MousePointer = 11 
Definitions.Show 
Definitions. WindowState = 0 


Screen.MousePointer = 0 
End Sub 


Sub CoolerName_Click () 
Screen.MousePointer = 11 
Cooler.Show 
Cooler. WindowsState = 0 
Screen.MousePointer = 0 

End Sub 


Sub CascadeName_Click () 
Screen.MousePointer = 11 
Cascade.Show 
Cascade. WindowState = 0 
Screen.MousePointer = 0 

End Sub 


Sub SinglePerformanceName_Click () 


Screen.MousePointer = 11 
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ScrollUpdate = False 


Tc = Val(Performance.LowTempText.Text) 
DeltaT Max = .5* z * (Tc ^2) 
TMax = DeltaTMax + Tc 


Performance.Tempscroll.Max = Tc * 10+ 1 
Performance.Tempscroll.Min = TMax * 10 - 1 


Performance. Tempscroll. Value = Performance.Tempscroll.Min 


Performance. TMinDisplay. Text = Format$(Tc, "##,##0.0") 
Performance. TMaxDisplay.Text = Format$(TMax, "##,##0.0") 


Th = Performance. Tempscroll. Value / 10 
Performance.HighTempDisplay. Text = Format$(Th, "##,##0.0") 
DeltaT = Th - Tc 

Performance.DeltaTDisplay.Text = Format$(DeltaT, "##,##0.0") 


Performance.Show 
Performance. WindowState = 0 


Screen.MousePointer = 0 
ScrollUpdate = True 

End Sub 

Sub MaterialName_Click () 
Screen.MousePointer = 11 


alpha_A = .00023 
alpha_B = .00021 


rho À = .001! 
rho B -» .001' 
kappa, A - .017 


kappa. B = .0145 


Aa, over, Àb 7 Sqr((kappa, B * rho. A) / (kappa, A * rho B)) 
Ab = .3848 

Aa= Aa_over_Ab * Ab 

dB = 2 * Sqr(Ab / pi) 

dA = 2 * Sqr(Aa / pi) 
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Par 
LoS Ea 
alpha = alpha_A + alpha_B 


z = (alpha ^ 2) / ((Sqr(kappa_A * rho_A) + Sqr(kappa_B * rho_B)) ^2) 
R 7 (rho. A * La / Aa) + (rho_B * Lb/ Ab) 
K = ((kappa_A * Aa) / La) + ((kappa_B * Ab) / Lb) 


Material.Optionl.Value = True 


Material.AlphaText.Text = Format$(alpha_A, "0.00000") 
Material.RhoText.Text = Format$(rho_A, "0.00000") 
Material.kText.Text = Format$(kappa_A, "0.00000") 


Material.FirstMaterial_a_Display.Text = Format$(Aa, "##0.0000") 
Material.FirstMaterial_d_Display.Text = Format$(dA, "##0.0000") 
Material.SecondMaterial_d_Display.Text = Format$(dB, "##0.0000") 
Material.FirstMaterial_L_Display.Text = Format$(La, "##0.0000") 
Material.RDisplay.Text = Format$(R, "0.000E+00") 
Material.KDisplay.Text = Format$(K, "0.000E+00") 
Material.zDisplay. Text = Format$(z, "0.000E+00") 


Materıal.Show 

Material. WindowState = 0 

Screen.MousePointer = 0 
End Sub 


Material.txt 
// Code for Material Properties Window 
Dim Aa, Ab, La, Lb As Double 


Sub Form_Load () 
Material.Left = (Screen. Width - Material. Width) / 2 
Material.Top 7 (Screen.Height - Material.Height) / 2 
End Sub 


Sub Materials Click () 
End Sub 


Sub Option! Click () 
If Optionl.Value = True Then 
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AlphaText.Text = Format$(alpha_A, "0.00000") 
RhoText.Text = Format$(rho_A, "0.00000") 
kText.Text = Format$(kappa_A, "0.00000") 
Else 
AlphaText.Text = Format$(alpha_B, "0.00000") 
RhoText.Text = Format$(rho_B, "0.00000") 
kText.Text = Format$(kappa_B, "0.00000") 
End If 
End Sub 


Sub Option2_Click () 

If Option]. Value = True Then 
AlphaText.Text = Format$(alpha_A, "0.00000") 
RhoText.Text = Format$(rho_A, "0.00000") 
kText.Text = Format$(kappa_A, "0.00000") 

Else 
AlphaText.Text = Format$(alpha_B, "0.00000") 
RhoText.Text = Format$(rho_B, "0.00000") 
kText.Text = Format$(kappa_B, "0.00000") 

End If 


End Sub 


Sub AlphaText_Change () 
If Option]. Value = True Then 
alpha_A = Val(AlphaText.Text) 
Else 
alpha_B = Val(AlphaText.Text) 
End If 
End Sub 


Sub RhoText_Change () 
If Option]. Value = True Then 
rho_A = Val(RhoText.Text) 
Else 
rho_B = Val(RhoText.Text) 
End If 
End Sub 


Sub kText Change () 
If Optionl.Value = True Then 
kappa_A = Val(kText.Text) 
Else 
kappa_B = Val(kText.Text) 
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End If 
End Sub 


Sub Commandl, Click () 
Main.OptimizeName.Enabled = True 
Material. Hide 

End Sub 


Sub Command2_Click () 
alpha_A = .00023 
alpha_B = .00021 


rho_A = .001! 
rho_B = .001! 
kappa_A =.017 


kappa, B = .0145 


Aa over Ab - Sqr((kappa, B * rho A) / (kappa, A * rho. B)) 
Ab = .3848 

Aa = Aa_over_Ab * Ab 

dB = 2 * Sqr(Ab / pi) 

dA = 2 * Sqr(Aa/ pi) 

La = .318 

Lb = La 

alpha = alpha_A + alpha_B 


z = (alpha ^ 2) / ((Sqr(kappa_A * rho_A) + Sqr(kappa_B * rho_B)) ^2) 
R = (rho_A * La / Aa) + (rho_B * Lb / Ab) 
K = ((kappa_A * Aa) / La) + ((kappa_B * Ab) / Lb) 


Option1.Value = True 


AlphaText.Text = Format$(alpha_A, "0.00000") 
RhoText.Text = Format$(rho_A, "0.00000") 
kText.Text = Format$(kappa_A, "0.00000") 


FirstMaterial_a_Display.Text = Format$(Aa, "##0.0000") 
SecondMaterial_a_Display.Text = Format$(Ab, "##0.0000") 
FirstMaterial_d_Display.Text = Format$(dA, "##0.0000") 
SecondMaterial_d_Display.Text = Format$(dB, "##0.0000") 
FirstMaterial_L_Display.Text = Format$(La, "##0.0000") 
SecondMaterial_L_Display.Text = Format$(Lb, "##0.0000") 
RDisplay.Text = Format$(R, "0.000E+00") 

KDisplay.Text = Format$(K, "0.000E+00") 

zDisplay.Text = Format$(z, "0.000E+00") 
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End Sub 


Sub Command3, Click () 
Aa, over, Ab 7 Sqr((kappa, B * rho, A) / (kappa, A * rho. B)) 
Ab = Val(SecondMaterial_a_Display.Text) 
Aa ^ Aa over Ab * Ab 
dB = 2 * Sqr(Ab / pi) 
dA = 2 * Sqr(Aa / pi) 
Lb = Val(SecondMaterial_L_Display.Text) 
La=Lb 
alpha = alpha_A + alpha_B 


z = (alpha ^ 2) / ((Sqr(kappa, A * rho A) + Sqr(kappa_B * rho_B)) ^ 2) 
R 7 (rho A * La/ Aa) * (rho B* Lb/ Ab) 
K 7 ((kappa, A * Aa) / La) * ((kappa B * Ab)/ Lb) 


FirstMaterial_a_Display. Text = Format$(Aa, "##0.0000") 
SecondMaterial_a_Display.Text = Format$(Ab, "##0.0000") 
FirstMaterial_d_Display.Text = Format$(dA, "##0.0000") 
SecondMaterial_d_Display.Text = Format$(dB, "##0.0000") 
FirstMaterial_L_Display. Text = Format$(La, "##0.0000") 
SecondMaterial_L_Display. Text = Format$(Lb, "##0.0000") 
RDisplay.Text = Format$(R, "0.000E+00") 

KDisplay.Text = Format$(K, "0.000E+00") 

zDisplay.Text = Format$(z, "0.000E+00") 


End Sub 


Coinput.txt 


Sub Command1_Click () 
Screen.MousePointer = 11 
Stage = 0 
Performance.Hide 
History.Hide 
Results.Hide 
Screen.MousePointer = 0 

End Sub 


Sub Form_Load () 


PerformanceCmd. Width = 2415 
Performance.Left = Screen. Width - Performance. Width 
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Performance.Top = Main.Height 
LowTempUpdate = True 
Cascading = False 


End Sub 
Sub Command2 Click () 


Screen.MousePointer = 11 
Command2.Caption = "Wait" 
History. WindowState — 0 
LowTempUpdate = True 
Cascading = True 


Stage = Stage + 1 


If Stage = 1 Then 
TotalPower = 0 
End If 


Tc = Val(LowTempText.Text) 
Th = Val(HighTempDisplay.Text) 


DeltaT = Th - Tc 
Th_plus_Tc = Th + Te 
gamma = Sqr(1 + (z / 2) * (Th, plus Tc)) 
Io = (alpha * DeltaT) / (R * (gamma - 1)) 
Vo = (alpha * DeltaT) + lo *R 
P -Io* Vo 
qT = Val(ReferLoadText.Text) 
qP = alpha * lo * Tc 
qJ =.5* (Io^2)*R 
qF = K * DeltaT 
qo = qP - (qJ + qF) 
No_of_Couples = qT / qo 
Total_Heat_rejected = qT + No_of_Couples * P 
Coeff_of_Performance = qo / P 
TotalPower = TotalPower + (P * No_of_Couples) 
If Stage = 1 Then 
FirstStage_qo = qo * No_of_Couples 
End If 
OverallCOP = FirstStage_qo / TotalPower 


TMaxDisplay.Text = Format$(DeltaTMax + Tc, "##,##0.0") 
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TMinDisplay.Text = Format$(Tc, "##,##0.0") 
DeltaTDisplay.Text = Format$(DeltaT, "##,##0.0") 


Results.CoeffText.Text = Format$(Coeff_of_Performance, "0.0000") 

Results. NumberCouplesText.Text = Format$(No_of_Couples, "#,4#0.00") 
Results. TotalHeatRejectedText.Text = Format$(Total_Heat_rejected, "#,4##0.0") 
Results.AperCoupleText.Text = Format$(Io, "#,##0.00") 

Results. VperCoupleT ext.Text = Format$( Vo, "#,4##0.0000") 
Results.WperCoupleText.Text = Format$(P, "#,4##0.0000") 

Results.qTText.Text = Format$(qT, "#,##0.0") 

Results.qPText.Text = Format$(qP, "#,##0.0000") 

Results.qJText. Text = Format$(qJ, "#,4##0.0000") 

Results.qF Text.Text = Format$(qF, "#,##0.0000") 

Results.qoText.Text = Format$(qo, "#,4##0.0000") 

Results.Caption = "Optimum Coefficient of Performance" + " Stage " + Str$(Stage) 
Results.Show 


Select Case Stage 
Case | 

History.qLoadStageOneDisplay.Text = Format$(qT, "#,##0.0") 

History.qRejectedStageOneDisplay.Text = Format$(Total_Heat_rejected, 
"# ##0.0") 

History.CouplesStageOneDisplay. Text = Format$(No_of_Couples, 
"#,##0.00") 

History. WperCoupleStageOneDisplay.Text = Format$(P, "#,##0.0000") 

History.TcStageOneDisplay.Text = Format$(Tc, "0.0000") 

History. ThStageOneDisplay.Text = Format$(Th, "0.0000") 

History.DeltaTStageOneDisplay.Text = Format$(DeltaT, "0.0000") 

History.StageCOPStageOneDisplay.Text = Format$(Coeff_of_Performance, 
"0.0000") 

History.OverallCOPStageOneDisplay.Text = Format$(OverallCOP, "0.0000") 

History. Width = 2250 

Case 2 

History.qLoadStageTwoDisplay.Text = Format$(qT, "#,##0.0") 

History.qRejectedStageT woDisplay.Text —^ Format$(Total Heat rejected, 
"#,##0.0") 

History.CouplesStageT woDisplay.Text ^ Format$(No, of. Couples, 
"#,##0.00") 

History. WperCoupleStageTwoDisplay.Text = Format$(P, "#,##0.0000") 

History.TcStageTwoDisplay.Text = Format$(Tc, 0.0000") 

History. ThStageTwoDisplay.Text = Format$(Th, "0.0000") 

History.DeltaTStageTwoDisplay. Text = Format$(DeltaT, "0.0000") 

History.StageCOPStageT woDisplay.Text ^ Format$(Coeff of Performance, 
"0.0000") 
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History.OverallCOPStageTwoDisplay. Text = Format$(OverallCOP, "0.0000") 
History. Width = 3360 
Case 3 
History.gLoadStageThreeDisplay.Text = Format$(qT, "#,##0.0") 
History.qRejectedStageThreeDisplay.Text ^ Format$(Total Heat rejected, 
"#,##0.0") 
History.CouplesStageThreeDisplay. Text = Format$(No_of_Couples, 
"#,##0.00") 
History. WperCoupleStageThreeDisplay. Text = Format$(P, "#,##0.0000") 
History.TcStageThreeDisplay.Text = Format$(Tc, "0.0000") 
History. ThStageThreeDisplay.Text = Format$(Th, "0.0000") 
History.DeltaTStageThreeDisplay.Text = Format$(DeltaT, "0.0000") 
History.StageCOPStageThreeDisplay.Text ^ Format$(Coeff of Performance, 
"0.0000") 
History.OverallCOPStageThreeDisplay.Text = Format$(OverallCOP, 
"0.0000") 
History. Width = 4425 
Case 4 
History.qLoadStageFourDisplay.Text = Format$(qT, "#,##0.0") 
History.qRejectedStageFourDisplay. Text = Format$(Total_Heat_rejected, 
"#,##0.0") 
History.CouplesStageFourDisplay.Text ^ Format$(No. of. Couples, 
"#,##0.00") 
History.WperCoupleStageFourDisplay. Text = Format$(P, "#,##0.0000") 
History. TcStageFourDisplay. Text = Format$(Tc, "0.0000") 
History. ThStageFourDisplay.Text = Format$(Th, "0.0000") 
History.DeltaTStageFourDisplay. Text = Format$(DeltaT, "0.0000") 
History.StageCOPStageFourDisplay.Text —^ Format$(Coeff of Performance, 
"0.0000") 
History.OverallCOPStageFourDisplay. Text = Format$(OverallCOP, "0.0000") 
History. Width = 5520 
Case 5 
History.qLoadStageFiveDisplay.Text = Format$(qT, "#,##0.0") 
History.qRejectedStageFiveDisplay.Text = Format$(Total_Heat_rejected, 
"#,##0.0") 
History.CouplesStageFiveDisplay.Text = Format$(No_of_Couples, 
"#,##0.00") 
History. WperCoupleStageFiveDisplay.Text = Format$(P, "#,##0.0000") 
History. TcStageFiveDisplay. Text = Format$(Tc, "0.0000") 
History. ThStageFiveDisplay.Text = Format$(Th, "0.0000") 
History.DeltaTStageFiveDisplay.Text = Format$(DeltaT, "0.0000") 
History.StageCOPStageFiveDisplay.Text = Format$(Coeff_of_Performance, 
"0.0000") 
History.OverallCOPStageFiveDisplay. Text = Format$(OverallCOP, "0.0000") 
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History.Width = 6585 
End Select 


History.Show 0 


ScrollUpdate = False 
LowTempText.Text = Format$(Th, "##,##0.0") 
ScrollUpdate = True 


ReferLoadText.Text = Format$(Total_Heat_rejected, "##0.0") 


LowTempUpdate = False 
Tempscroll. Value = (Tc + DeltaT) * 10 
LowTempUpdate = True 


Command2.Caption = "Cascade" 
Screen.MousePointer = 0 


Cascading = False 
End Sub 


Sub TempScroll_Change () 
If ScrollUpdate = True Then 
Tc = Val(LowTempText.Text) 
Th = Tempscroll. Value / 10 
DeltaT = Th - Tc 
HighTempDisplay.Text = Format$(Th, "##,#00.0") 
DeltaTDisplay.Text = Format$(DeltaT, "##,#00.0") 
End If 
End Sub 


Sub LowTempText Change () 

If LowTempUpdate = True Then 
Tc = Val(LowTempText.Text) 
DeltaTMax = .5 *z* (Tc ^ 2) 
1 Max 7 Deita1Max + tc 


Tempscroll.Max = Tc * 10 + 1 
Tempscroll.Min = TMax * 10 - 1 
If Stage = 0 Then 

Tempscroll. Value = Tempscroll.Min 
End If 


Th = Tempscroll. Value / 10 


53 


If Cascading = False Then 
DeltaT = Th - Tc 
End If 


TMinDisplay.Text = Format$(Tc, "##,##0.0") 
TMaxDisplay.Text = Format$(TMax, "##,##0.0") 
DeltaTDisplay. Text = Format$(DeltaT, "##,##0.0") 
End If 
End Sub 


Sub PerformanceCmd_Click () 


Screen.MousePointer = 11 
Main.ShowGraphName.Enabled = True 
PerformanceCmd.Caption — "Wait" 
Command2. Visible = -1 

History.Hide 


Stage = 0 


DeltaT = Th - Tc 

Th plus Tc = Th + Tc 

gamma = Sar(1 +(z/2) ” (Th plus. Tc)) 
Io = (alpha * DeltaT) / (R * (gamma - 1)) 
Vo = (alpha * DeltaT) + Io *R 

P-Io * Vo 

qT = Val(ReferLoadText.Text) 

qP - alpha * Io * Tc 

a) = Sao DER 

qF=K * DeltaT 

qo 7 qP - (qJ * qF) 

No of Couples = qT / qo 
Total_Heat_rejected = qT + No_of_Couples * P 
Coeff_of_Performance = qo / P 


TMaxDisplay.Text = Format$(DeltaTMax + Tc, "##,##0.0") 
TMinDisplay.Text = Format$(Tc, "##,##0.0") 
DeltaTDisplay.Text = Format$(DeltaT, "##,##0.0") 


Results.CoeffText.Text = Format$(Coeff_of_Performance, "0.0000") 
Results.NumberCouplesText. Text = Format$(No_of_Couples, "#,##0.00") 
Results. TotalHeatRejectedText.Text = Format$(Total_Heat_rejected, "#,##0.0") 
Results. AperCoupleText.Text = Format$(Io, "#,##0.00") 
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Results. VperCoupleText.Text = Format$(Vo, "#,##0.0000") 
Results. WperCoupleText.Text = Format$(P, "#,##0.00") 
Results.qTText.Text = Format$(qT, "#,##0.0") 
Results.qPText.Text = Format$(qP, "#,4##0.0000") 
Results.qJText.Text = Format$(qJ, "#,##0.0000") 
Results.qFText.Text = Format$(qF, "#,##0.0000") 
Results.goText.Text = Format$(go, "#,##0.0000") 
Results.Caption = "Optimum Coefficient of Performance" 
Results.Show 
Results.PowerConsumptionDisplay.Text = Format$((P * No_of_Couples), 
"# #H# ##0.0") 
PerformanceCmd. Width = 975 
PerformanceCmd.Caption - "Restart" 
Screen.MousePointer = 0 
End Sub 


Cooler.txt 
// Code for Optimum Coeficient of Performance Window 
DefDbl A-Z 


Global Const True = -1 
Global Const False = 0 
Global Const pi = 355/113 
Global DeltaT As Double 
Global DeltaTMax As Double 


Global alpha_A, alpha_B, rho_A, rho_B, kappa_A, kappa_B As Double 
Global alpha As Double 

Global z, K, R As Double 

Global Stage As Integer 

Global ScrollUpdate As Integer 
Global LowTempUpdate As Integer 
Global Cascading As Integer 

Global Total Heat rejected As Double 
Global Tc As Double 

Global Th As Double 

Global TotalPower As Double 

Global FirstStage qo As Double 


Cascade.txt 

// Code for Cascaded Stages Window 

Sub PrintName Click () 
Cascade.PrintForm 
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End Sub 


Sub ExitName Click () 
Cascade.Hide 
End Sub 


Define.txt 

// Code to call Definitions Window 

Sub PrintName Click () 
Definitions.PrintForm 

End Sub 


Sub ExitName Click () 
Definitions.Hide 
End Sub 


Diagram.txt 

// Code to call Diagrams 

Sub PrintName, Click () 
Cooler.PrintForm 

End Sub 


Sub ExitName, Click () 
Cooler.Hide 
End Sub 


Ending.txt 
// Code to end design session, closing all open windows 
Sub Commandl. Click () 
End 
End Sub 


Sub Command2, Click () 
Ending.Hide 
End Sub 
Sub Form, Click () 
Main.Left = (Screen. Width - Main. Width) / 2 
Main.Top = (Screen.Height - Main.Height) / 2 
End Sub 


Sub Form_Load () 
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Ending.Left = (Screen. Width - Ending. Width) / 2 
Ending.Top = (Screen.Height - Ending.Height) / 2 


End Sub 


Graph.txt 

// Code for graping routines 

Dim Shown As Integer 

Dim Samples As Integer 

Dim IMax As Double 

Dim IMin As Double 

Dim Y Max As Double 

Dim YMin As Double 

Dim OldY Max As Double 

Dim OldIMax As Double 

Dim OldIMin As Double 

Dim OldMaxCount As Double 
Dim OldMinCount As Double 
Dim IatOldMinCount As Double 
Dim latOldMaxCount As Double 
Dim PerformancePlot As Integer 


Dim SuperimposePlots As Integer 


Sub Graphlt () 


ReDim I(Samples) 
ReDim Y(Samples) 


IMin = ((2 * alpha * Tc) - Sqr((2 * alpha * Tc) ^2 -(8 * R * K * DeltaT))) / (2 * R) 
[Max alpha " Tc) + Sqr((2 * alpha * Tc) "2 -(8 RÆK * DeltaT))) 4(2 * R) 


YMin = 0 
YMax = 0 
MinCount = 0 
MaxCount = 0 


Minimum = False 


StepInterval = (IMax - IMin) / (Samples - 1) 


If PerformancePlot = True Then 
Y AxisLabel Text = "COP" 
Graph.Caption = "Performance" 


Else 
Y AxisLabel. Text = "q" 


Graph.Caption = "Heat Pumping" 


B 


End If 


For J= 1 To Samples 
I(J) = IMin + StepInterval * (J - 1) 


If PerformancePlot = True Then 
Y (J) = (Calpha * 11) * TO) -(S5*I0)*2*R)-(K* DIED ap 
DeltaT) + (1(J) ” R))) 
Else 
Y(J) = (alpha * Tc * I(J)) - (.5 + 10) “2 £ R) - (K £ DeltaT) 
End If 


If Y(J) > 0 And Minimum = False Then 
YMin = Y(J) 
MinCount = J 
Minimum = True 
Elself Y(J) > YMax Then 
Y Max = Y(J) 
| at YMax - I(J) 
Elself Y(J) > 0 Then 
MaxCount = J 
End If 


Next J 


If SuperimposePlots = True Then 
ResultsPicture.ScaleMode = 0 
ResultsPicture.ScaleTop = OldY Max 
ResultsPicture.ScaleHeight = -OldY Max 
ResultsPicture.ScaleLeft = latOldMinCount 
ResultsPicture.ScaleWidth 7 IatOldMaxCount - IatOldMinCount 
YMaxText.Text = Format$(OldYMax, "0.00 e+00") 
IMinText.Text = Format$(OldIMin, "0.00 e+00") 
[MaxText.Text = Format$(OldIMax, "0.00 e+00") 

Else 
ResultsPicture.Cls 
ResultsPicture.ScaleMode = 0 
ResultsPicture.ScaleTop = YMax 
ResultsPicture.ScaleHeight = -Y Max 
ResultsPicture.ScaleLeft = [(MinCount) 
ResultsPicture.ScaleWidth = I(MaxCount) - I(MinCount) 
Y MaxText.Text 7 Format$( Y Max, "0.00 e*00") 
IMinText.Text = Format$(IMin, "0.00 e+00") 
IMaxText. Text = Format$(IMax, "0.00 e+00") 
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OldYMax = YMax 

OldIMax = IMax 

OldIMin = IMin 

OldMaxCount = MaxCount 

IatOldMinCount 7 I( MinCount) 

IatOldMaxCount = I(MaxCount) 
End If 


For J = 1 To Samples 
ResultsPicture.Line (I(J - 1), Y(J - 1))-(I(J), Y(J)) 
Next J 


I at. YMaxDisplay.Text ^ Format$(I at YMax, "0.00 e+00") 
DeltaTDisplay.Text = Format$(DeltaT, "0.00") 


End Sub 


Sub ExitName_Click () 
Graph.Hide 
End Sub 


Sub PrintName, Click () 
Graph.PrintForm 
End Sub 


Sub PerformanceName Click () 
PerformancePlot = True 
PerformanceName.Checked = True 
PumpName.Checked = False 

End Sub 


Sub PumpName, Click () 
PerformancePlot = False 
PerformanceName.Checked = False 
PumpName.Checked = True 

End Sub 


Sub Form_Load () 
PerformancePlot = True 
Samples = 200 

End Sub 
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Sub DrawName Click () 
Screen.MousePointer 7 11 
Graphlt 
SuperName.Enabled = True 
Screen.MousePointer = 0 

End Sub 


Sub al00Name Click () 
al00Name.Checked = True 
a200Name.Checked = False 
aS00Name.Checked = False 
al000Name.Checked = False 
Samples = 100 

End Sub 


Sub a200Name_Click () 
al00Name.Checked = False 
a200Name.Checked = True 
aS00Name.Checked = False 
al000Name.Checked = False 
Samples = 200 


End Sub 


Sub a500Name Click () 
al00Name.Checked = False 
a200Name.Checked = False 
a500Name.Checked = True 
al000Name.Checked = False 
Samples = 500 


End Sub 


Sub al000Name_Click () 
a100Name.Checked = False 
a200Name.Checked = False 
a500Name.Checked = False 
al000Name.Checked 7» True 
Samples = 1000 


End Sub 


Sub SuperName, Click () 
If SuperName.Checked = False Then 


60 


SuperimposePlots = True 
SuperName.Checked = True 
Else 
SuperimposePlots = False 
SuperName.Checked = False 
End If 
End Sub 


History.txt 
// Cascaded Stages Window menu bar 
Sub Form_Load () 
History.Left = 0 
History.Top = Main.Height 
End Sub 


Sub PrintName, Click () 
History.PrintForm 
End Sub 


Intro.txt 
// Introduction graphic, program startup 
Sub Timerl, Timer () 
Screen.MousePointer = 11 
Load Results 
Load Performance 
Load Material 
Load Definitions 
Load TempCalc 
Load Graph 
Load Cascade 
Main.Show 
Load Cooler 
Load Ending 
Load History 
Intro. Hide 
Tımerl.Enabled = 0 
Screen.MousePointer = 0 
End Sub 


Sub Form_Load () 


Screen.MousePointer = 11 
End Sub 
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Pumping.txt 
// Code for Maximum Heat Pumped Window 
Sub Commandl. Click () 
Screen.MousePointer 7 11 
Stage = 0 
Pumping.Hide 
History.Hide 
Results. Hide 
Screen.MousePointer = 0 
End Sub 


Sub Form_Load () 


PerformanceCmd. Width = 2415 

Pumping.Left = Screen.Width - Pumping. Width 
Pumping.Top = Main.Height + 450 
LowTempUpdate = True 

Cascading = False 


End Sub 
Sub Command2 Click () 


Screen.MousePointer = 11 
Command2.Caption = "Wait" 
History. WindowState = 0 
LowTempUpdate = True 
Cascading = True 


Stage = Stage + 1 


If Stage = 1 Then 
TotalPower = 0 
End If 


Tc = Val(LowTempText.Text) 
Th = Val(HighTempDisplay.Text) 


DeltaT = Th - Tc 

Delta MAS ic) 
Th_plus_Tc = Th + Tc 

gamma = Sqr(1 + (z/ 2) * (Th plus Tc)) 
lo = (alpha * Tc) /R 
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Vo = alpha * Th 
P = [o * Vo 
qT = Val(ReferLoadText.Text) 
qP = alpha * Io * Tc 
11 068 6)” ከ 
qF » K * DeltaT 
qo 7.5 *((alpha^2 * Tc^2)/ R) - gF 
No. of Couples = qT / qo 
Total_Heat_rejected = qT + No_of_Couples * P 
Coeff of Performance 7 qo / P 
TotalPower = TotalPower * (P * No of Couples) 
If Stage = 1 Then 
FirstStage qo ^ qo * No of Couples 
End If 
OverallCOP 7 FirstStage qo / TotalPower 


TMaxDisplay.Text = Format$(DeltaTMax + Tc, "##,##0.0") 
TMinDisplay.Text = Format$(Tc, "##,##0.0") 
DeltaTDisplay. Text = Format$(DeltaT, "##,##0.0") 


Results.CoeffText.Text ^ Format$(Coeff of Performance, "0.0000") 
Results. NumberCouplesText.Text = Format$(No_of_Couples, "#,##0.00") 
Results. TotalHeatRejectedText.Text = Format$(Total_Heat_rejected, "#,4##0.0") 
Results. AperCoupleText.Text = Format$(Io, "#,4##0.00") 

Results. VperCoupleText. Text = Format$( Vo, "#,##0.0000") 

Results. WperCoupleText. Text = Format$(P, "#,##0.0000") 
Results.qTText.Text = Format$(qT, "#,##0.0") 

Results.qPText.Text = Format$(qP, "#,4##0.0000") 

Results.qJText.Text = Format$(qJ, "#,##0.0000") 

Results.qFText.Text = Format$(qF, "#,##0.0000") 

Results.goText.Text = Format$(go, "#,##0.0000") 

Results.Caption = "Maximum Heat Pumping" + " Stage " + Str$(Stage) 
Results.Show 


Select Case Stage 
Case | 

History.qLoadStageOneDisplay.Text = Format$(qT, "#,##0.0") 

History.qRejectedStageOneDisplay.Text = Format$(Total_Heat_rejected, 
"#,##0.0") 

History.CouplesStageOneDisplay. Text = Format$(No_of_Couples, 
"#,##0.00") 

History. WperCoupleStageOneDisplay.Text = Format$(P, "#,##0.0000") 

History.TcStageOneDisplay.Text = Format$(Tc, "0.0000") 

History. ThStageOneDisplay.Text = Format$(Th, "0.0000") 
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History.DeltaTStageOneDisplay. Text = Format$(DeltaT, "0.0000") 
History.StageCOPStageOneDisplay.Text ^ Format$(Coeff of Performance, 
"0.0000") 
History.OverallCOPStageOneDisplay.Text = Format$(OverallCOP, "0.0000") 
History. Width = 2250 
Case 2 
History.qLoadStageT woDisplay.Text = Format$(qT, "#,##0.0") 
History.qRejectedStageT woDisplay.Text = Format$(Total_Heat_rejected, 
"#,##0.0") 
History.CouplesStageT woDisplay.Text ^ Format$(No, of Couples, 
"#,##0.00") 
History. WperCoupleStageTwoDisplay.Text = Format$(P, "#,##0.0000") 
History.TcStageT woDisplay.Text = Format$(Tc, "0.0000") 
History. ThStageT woDisplay.Text = Format$(Th, "0.0000") 
History.DeltaTStageT woDisplay.Text ^ Format$(DeltaT, "0.0000") 
History.StageCOPStageT woDisplay.Text — Format$(Coeff of Performance, 
"0.0000") 
History.OverallCOPStageTwoDisplay.Text = Format$(OverallCOP, "0.0000") 
History. Width = 3360 
Case 3 
History.qLoadStageThreeDisplay.Text = Format$(qT, "#,##0.0") 
History.qRejectedStageThreeDisplay.Text — Format$(Total Heat rejected, 
"f 820.0") 
History.CouplesStageThreeDisplay.Text ^ Format$(No. of Couples, 
"£ 880.00") 
History.WperCoupleStageThreeDisplay.Text = Format$(P, "#,##0.0000") 
History.TcStageThreeDisplay.Text = Format$(Tc, "0.0000") 
History. ThStageThreeDisplay.Text = Format$(Th, "0.0000") 
History.DeltaTStageThreeDisplay.Text = Format$(DeltaT, "0.0000") 
History.StageCOPStageThreeDisplay.Text = Format$(Coeff_of_Performance, 
"0.0000") 
History.OverallCOPStageThreeDisplay.Text = Format$(OverallCOP, 
"0.0000") 
History. Width = 4425 
Case 4 
History.qLoadStageFourDisplay.Text = Format$(qT, "#,##0.0") 
History.qRejectedStageFourDisplay.Text = Format$(Total_Heat_rejected, 
"#,##0.0") 
History.CouplesStageFourDisplay.Text = Format$(No_of_Couples, 
"#,##0.00") 
History.WperCoupleStageFourDisplay.Text = Format$(P, "#,##0.0000") 
History.TcStageFourDisplay.Text = Format$(Tc, "0.0000") 
History.ThStageFourDisplay.Text = Format$(Th, "0.0000") 
History.DeltaTStageFourDisplay. Text = Format$(DeltaT, "0.0000") 
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History.StageCOPStageFourDisplay.Text ^ Format$(Coeff of Performance, 
"0.0000") 

History.OverallCOPStageFourDisplay.Text = Format$(OverallCOP, "0.0000") 

History. Width = 5520 

Case 5 

History.qLoadStageFiveDisplay.Text = Format$(qT, "#,##0.0") 

History.qRejectedStageFiveDisplay. Text = Format$(Total_Heat_rejected, 
"ዴ##0.0") 

History.CouplesStageFiveDisplay.Text ^ Format$(No of Couples, 
"#,##0.00") 

History. WperCoupleStageFiveDisplay.Text = Format$(P, "#,##0.0000") 

History. TcStageFiveDisplay. Text = Format$(Tc, "0.0000") 

History. ThStageFiveDisplay. Text = Format$(Th, "0.0000") 

History.DeltaTStageFiveDisplay.Text = Format$(DeltaT, "0.0000") 

History.StageCOPStageFiveDisplay. Text = Format$(Coeff_of_Performance, 


"0.0000") 
History.OverallCOPStageFiveDisplay. Text = Format$(OverallCOP, "0.0000") 
History. Width = 6585 
End Select 
History.Show 0 


ScrollUpdate = False 
LowTempText.Text = Format$(Th, "##,##0.0") 
ScrollUpdate = True 


ReferLoadText.Text = Format$(Total_Heat_rejected, "##0.00") 


LowTempUpdate = False 
Tempscroll. Value = (Tc + DeltaT) * 10 
LowTemp Update = True 


Command2.Caption = "Cascade" 
Screen.MousePointer = 0 


Cascading = False 
End Sub 


Sub TempScroll_Change () 
If ScrollUpdate = True Then 
Te = Val(LowTempText.Text) 
Th = Tempscroll. Value / 10 
DeltaT = Th - Tc 
HighTempDisplay.Text = Format$(Th, "##,#00.0") 
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DeltaTDisplay.Text = Format$(DeltaT, "##,#00.0") 
End If 
End Sub 


Sub LowTempText_Change () 

If LowTempUpdate = True Then 
16 = Val(LowTempText.Text) 
DeltaTMax = .5 * z * (Tc ^2) 
TMax = DeltaTMax + Tc 


Tempscroll.Max = Tc * 10 + 1 
Tempscroll. Min = TMax * 10 - 1 
If Stage = 0 Then 

Tempscroll. Value = Tempscroll.Min 
End If 


Th = Tempscroll. Value / 10 


If Cascading = False Then 
DeltaT = Th - Tc 
End If 


TMinDisplay.Text = Format$(Tc, "##,##0.0") 
TMaxDisplay.Text = Format$(TMax, "##,##0.0") 
DeltaTDisplay.Text = Format$(DeltaT, "##,##0.0") 
End If 
End Sub 


Sub PerformanceCmd Click () 


Screen. MousePointer = 11 
Main.ShowGraphName.Enabled = True 
PerformanceCmd.Caption = "Wait" 
Command2. Visible = -1 

History.Hide 


Stage = 0 
DeltaT = Th - Tc 
Th plus Tc » Th * Tc 
gamma = Sqr(1 + (z/ 2) * (Th plus Tc)) 


Io 7 (alpha * Tc) / R 
Vo = alpha * Th 
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P = Io * Vo 

qT = Val(ReferLoadText.Text) 

qP = alpha * Io * Tc 

gq) = Selon 2) ~* RK 

qF = K * DeltaT 

qo =.5 * ((alpha ^ 2 * Tc ^ 2)/ R) - qF 

No. of Couples = gT / qo 

Total Heat rejected = qT + No_of_Couples * P 
Coeff_of_Performance = qo / P 


TMaxDisplay.Text = Format$(DeltaTMax + Tc, "##,##0.0") 
TMinDisplay.Text = Format$(Tc, "##,##0.0") 
DeltaTDisplay.Text = Format$(DeltaT, "##,##0.0") 


Results.CoeffText. Text ^ Format$(Coeff of Performance, "0.0000") 

Results. NumberCouplesText.Text = Format$(No_of_Couples, "#,##0.00") 

Results. TotalHeatRejectedText.Text = Format$(Total_Heat_rejected, "#,##0.0") 

Results. AperCoupleText.Text = Format$(lo, "#,##0.00") 

Results. VperCoupleText.Text = Format$(Vo, "#,##0.0000") 

Results. WperCoupleText.Text = Format$(P, "#,4##0.00") 

Results.qT Text. Text = Format$(qT, "#,##0.0") 

Results.qPText.Text = Format$(qP, "#,##0.0000") 

Results.qJText.Text = Format$(qJ, "#,##0.0000") 

Results.qFText. Text = Format$(qF, "#,##0.0000") 

Results.goText. Text = Format$(go, "#,##0.0000") 

Results.Caption = "Maximum Heat Pumping" 

Results.Show 

Results.PowerConsumptionDisplay.Text ^ Format$((P * No. of Couples), 
".]ጩ,#60.00") 

PerformanceCmd.Width = 975 

PerformanceCmd.Caption = "Restart" 

Screen.MousePointer = 0 
End Sub 


Results.txt 

// Code for Results Window 

Sub Form Click () 
Results.Top = (Screen.Height - Results. Height) 
Results.Left ^ (Screen. Width - Results. Width) / 2 


End Sub 
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Sub Form Load () 
Results. Top = (Screen.Height - Results. Height) 
Results.Left = (Screen. Width - Results. Width) / 2 
End Sub 


Sub PrintName_Click () 
Results. PrintForm 
End Sub 


Tempcalc.txt 

// Code for Temperature Conversion Calculator 

Sub Command2 Click () 
ConvertDisplayText.Text — 
ConvertDisplayText.SetFocus 

End Sub 
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Sub ConvertCmd Click () 
Const True = -1 
Const False = 0 
Temperature = Val(ConvertDisplayText.Text) 
If FarenheitOutOption. Value = True Then 
If RankineInOption. Value = True Then 
Temperature = Temperature - 459.67 
Elself FarenheitInOption. Value = True Then 
Temperature = Temperature 
Elself CelsiusInOption. Value = True Then 
Temperature = (Temperature * 9 / 5) + 32 
ElseIf KelvinInOption. Value = True Then 
Temperature = ((Temperature - 273,15) * 9/5) + 32 
End If 


Elself RankineOutOption. Value = True Then 


If RankineInOption. Value = True Then 
Temperature = Temperature 
Elself FarenheitInOption. Value = True Then 
Temperature = Temperature + 459.67 
Elself CelsiusInOption. Value = True Then 
Temperature = ((Temperature) * 9 / 5) + 32 + 459.67 
ElseIf KelvinInOption. Value = True Then 
Temperature = ((Temperature - 273.15) * 9/5) + 32 + 459.67 
End If 
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Elself CelsiusOutOption. Value = True Then 


If RankineInOption. Value = True Then 
Temperature = (Temperature - 32 - 459.67) * 5/9 

Elself FarenheitInOption. Value = True Then 
Temperature — (Temperature - 32) * 5/9 

Elself CelsiusInOption. Value = True Then 
Temperature = Temperature 

Elself KelvinInOption. Value = True Then 
Temperature = Temperature - 273.15 

End If 


Elself KelvinOutOption. Value = True Then 


If RankineInOption. Value = True Then 
Temperature = (Temperature - 32 - 459.67) * 5/9 + 273.15 
Elself FarenheitInOption. Value = True Then 
Temperature — (Temperature - 32) * 5/9 + 273.15 
Elself CelsiusInOption. Value = True Then 
Temperature = Temperature + 273.15 
ElseIf KelvinInOption. Value = True Then 
Temperature = Temperature 
End If 


End If 

ConvertDisplayText.Text = Format$(Temperature, "##,##0.00") 
ConvertDisplayText.SetFocus 

End Sub 


Sub Form Click () 
ConvertDisplayText.SetFocus 
End Sub 


Sub Form, Load () 
TempCalc.Left — (Screen. Width - TempCalc.Width) / 2 
TempCalc.Top = (Screen.Height - TempCalc.Height) / 2 
End Sub 


Sub ExitName_Click () 


TempCalc.Hide 
End Sub 
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